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1.0 INTRODUCTION 


This report summarizes the results of experiments in Curie depth esti- 
mation from long wavelength magnetic anomalies using methodologies previ- 
ously developed in connection with satellite data, both Pogo-series (Mayhew 
et al , 1980; Mayhew, 1982a, b) and Magsat (Mayhew, 1982c). The heart of the 
work is equivalent-layer-type magnetization models derived by inversion of 
high-elevation, long wavelength magnetic anomaly data. The methodology is 
described in detail in the above references. The magnetizat''on model for 
the United States used in Part 1 of the present work is givm in Mayhew 
(1982c), and in Figure 1 of this report. 

The central goal of the analysis is to find a magnetization distribu- 
tion in a thin equivalent layer at the Eartii's surface having maximum 
detail while retaining physical significance, and giving rise to a synthe- 
tic anomaly field which makes a best fit to the observed field in a least 
squares sense. Multiplying the magnetization distribution by the layer 
thickness gives the distribution of the vertical integral of magnetization 
in the magnetic crust to within an indeterminant ambiguity in level 
(Mayhew, 1982b). Such a model can be transformed to a more physically 
meaningful model given some indeoendent constraints, as described below ard 
in Mayhew (1982b). The apparent magnetization contrast {m) in the equiva- 
lent layer is approximated using an array of dipoles distributed in aqual 
area at the Earth's surface. The dipoles are pointed in the direction of 
the main magnetic, field, which carries the implicit assumption that crustal 
magnetization is dominantly induced or viscous. A key element of the 
analysis is the determination of the closest possible dipole spacing giving 
a "stable" inversion to a solution having physical significance. This is 
accomplished by plotting the standard deviation of the solution parameters 
(i.e, the magnetization values associated with the dipoles) against their 
spatial separation for a series of solutions. At separations closer than 
some critical one, the magnetization values become large (positive and 
negative) and no longer contour systematically. For Pogo data, the resolu- 
tion limit is at about 300 km separation. This is to be compared with a 
similar result based on Magsat data which implies a resolution limit 
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approaching 200 km for this lower data. Mayhew (1982a) shows how to find a 
magnetization distribution on a grid twice as fine as that at the resolu- 
tion limit, and this technique has been used In producing the am model 
shown In Figure 1, 

One of the most straight forward approaches to Interpretation Is to 
attempt to convert the equivalent layer am models to models of thickness 
variation In a layer of constant magnetization. This Is applicable, for 
example, to the case In which magnetic field variations reflect undulations 
of the Curie Isotherm within the crust. The methodology Is described In 
Mayhew (1982b). The vertical integral of magnetization implied by the am 
model is H(M+am), where H is the equivalent layer thickness (arbitrarily 
taken to be 40 km), and M is the level ambiguity in the solution, assumed 
to be constant. For the case of thickness variation h In a layer of con- 
stant magnetization u , the vertical integral of magnetization is hu. 
These two cases are indistinguishable in the anomaly field at satellite 
elevation, so equating the two gives a basic equation 


H(M+Am) = yh. (1) 

For the Curie depth problem, if its depth can be estimated at at least two 
places (by use of thermal models or spectral estimates) the unknown M and 
M can be estimated, and h determined for the whole region from the im 
distribution. The implied Curie depth configuration can then be used to 
constrain a regional crustal geothermal model. Results of this type are 
given In Section 2. 

The resolution obtainable in such Curie depth models is limited by the 
data elevation. Mayhew (1982b) suggested that much higher resolution 
models could be obtained using upward continued aeromagnetic data in an 
analogous manner, where data coverage over a sufficiently large region is 
available. Part 2 of the present work is a test of this idea; results are 
described in Section 3.0. 
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To summarize the principal elements of the statement of work for this 
project, they are 1) to derive two>dimensiona1 Curie depth models for sev- 
eral areas of the western U.S. using the equivalent-layer magnetization 
model derived fran Magsat data, use the results to constrain finite-element 
thermal nwdels, and compare theoretical heat flov* predicted by the models 
with that inferred from direct measurements, and 2) to investigate the 
applicability of the Curie depth modeling methodology developed for satel- 
lite data to upward continued aircraft data. 
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2.0 PART 1; TWO-DIMENSIONAL THERMAL MODELS 


The production of regional heat flow maps Is problematical for a 
number of reasons. Heat flow data sets are Inherently "noisy" because of 
various transient, ground water, topographic, and other effects. Measure- 
ments are irregularly distributed and of variable quality. Conductive heat 
flow is difficult to separate from convective effects due to local move- 
ments of hydrothermal fl <ds or magmas in active areas. Curie isotherm 
undulations should be reflected in regional heat flow variations. To the 
extent that Am anomalies are due to Curie depth variations, they can be 
used along with heat flow measurements to guide the construction of region- 
al heat flow maps. Such an approach was used by Mayhew (1982c) to make a 
regional heat flow map for the western U.S. (Figure 2), which was shown to 
compare favorably with similar maps made by others based principally on 
borehole measurements (Blackwell, 1979) and on silica geothermometry 
(Swanberg et al, in press). The first part of the present work was an 
investigation of whether Am maps can be used in a more quantitative 
fashion to accurately predict heat flow and Curie depth. 

2.1 Approach 

Four sectiori lines were chosen for two-dimensional analysis (Figure 
1). These lines c"oss the boundaries of regional heat flow provinces in- 
ferred from measurements, characteristics of which are summarized briefly 
as follows. The whole of the Basin and Range is a heat flow high except 
for a low area in southern Nevada which may or may not reflect crustal heat 
conduction (Lachenbruch and Sass, 1977). The Basin and Range high extends 
around the southern margin of the Colorado Plateau (in a poorly-defined 
fashion) to tie Rio Grande Rift, clearly a linear zone of high heat flow, 
although measured . values scatter greatly. The Colorado Plateau has lower 
heat flow than surrounding regions, although generally higher than 
"normal". The axis of the Sierra Nevada has low heat flow. These gross 
characteristics show up in regional heat flow Tiaps, and are reflected in 
the Am map (Figure 1). Low Am values apparently reflect shallow Curie 
depths on a regional scale, and thus high conductive crustal heat flow. 



Business and Technological Systems, Inc 


Am variations along the four section lines were transformed to Curie 
Isotherm configuration in the following man'ier. The basic equation Is (1), 
from which variations In magnetic crustal thickness h are computed from 
variations in Am . Estimates are needed for the parameters M and u In 
equation (1), These are gotten by fitting a lint to a plot of estimates of 
range of h derived from spectral depth estimates, depth to 550®C (the 
assumed Curie temperature) given by thermal models constrained by surface 
heat flow measurements, or. In the case of the Sierra Nevada, seismic esti- 
mates of depth to mantle, against range of Am for the corresponding 
areas. The Moho is assumed to be the magnetic bottom for the Sierras, 
since the very low heat flow values measured along the axis of the belt 
imply depths to 550° well into the mantle, which is assumed to be non- 
magnetic (Wasilewski et al , 1979). The plot is given as Figure 3. The 
ranges of values are rather large, but give an overall trend in the right 
direction. Tvo lines are run through the data. One Is an eyeball best 
fit, and implies M and u of 0.72 and 1.11 A/m, respectively. The other 
implies M =0.95, u = 2.00 A/m, and was used mainly to give a comparison, 
although the values are similar to those inferred from a preliminary study 
of the couthern Rio Grande Rift (Mayhew, 1982b), Both pairs of parameters 
were used in the modeling. 

The Curie surface configurations (assumed to be isothermal and equal 
to 550° on all sections but one) were used to constrain a finite-element 
thermal model. This software was further developed from that used previ- 
ously (Mayhew, 1982b). Numerous runs were made. Most involved a simple 
one-layer model with surface heat production decreasing decrementally with 
depth from a surface value Aq, and constant thermal conductivity K . 
Values for Aq and K of 5 HPU and 6 (cgs), respectively, were used as 
starter values, and are those suggested by Diment et al (1975) as typical 
for the crust, although measured values on crystalline rocks have a great 
range. Other values were tried as well, as were two-layer models. It was 
found that the overall results are not profoundly affected by the choice of 
reasonable thermal parameters. 
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2.2 Results 


Figure 4 (Section 1) Is a section extending from the Pacific coast to 
the Great Plains, om Is the apparent magnetization contrast along the 
section from Figure 1. Also shown are the Curie depth surfaces (he) for 
the two pairs of (M,u) mentioned above and the computed heat flow profile 
associated with each thermal model. The circles represent heat flow meas- 
urements within about 100 km of the section. While they are characteris- 
tically scattered, they are In good general agreement with predicted heat 
flow, especially model A. The results seem to confirm the existence of 
zones of elevated heat flow In the Basin and Range adjacent to the Wasatch 
Front and along the northern extension of the Rio Grande Rift. 

Figure 5 (Section 2) crosses from the southern Basin and Range to the 
Colorado Plateau. Both models predict about the same high heat flow (3 
HFU) south of the Colorado Plateau. Model 2B predicts a heat flow level 
quite similar to that associated with the Plateau (about 1.5 HFU). The 
Curie surfaces are quite different; that for model A is at Moho depth 
beneath the Plateau, but well up into ti.e crust for model B. 

Section 3 (Figure 6) crosses tne central Rio Grande Rift. Model A was 
computed assuming one layer with Aq and K equal to 5.0 and 6.0, re- 
spectively. For comparison, a two-layer model assuming K = 6.5 in an 
upper layer 20 km thick and a lower layer with K = 5.0, and with heat pro- 
duction Aq equal to 5.3 in the upper layer and constant in the lower 
layer and equal to the decremental value at 20 km, (i.e. A 2 = Aq exp 
(-20/0), was computed. D, the decremental constant, is taken to be 10 km, 
the value generally assumed for the western U.S. The other parameters are 
like those favored by Blackwell (1971). The computed heat flow profile is 
very close to model Model B also assumed one layer with Ag = 5,0 and 
K * 6.0, but different (M,u). The Curie surface for this model is shown, 
Wasilewski (personal communication) has measured low Curie points on lower 
crustal xenoliths from the southern Rio Grande Rift (about 300“C) and has 
speculated that this is a general characteristic of the rift zone, and the 
effect of reducing conduction in a hot, dry environment (\t deep crustal 
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levels. More “normal" values (around 550*C) are measured for the Colorado 
Plateau. Model C is for a hypothetical case In which the Curie surface Is 
not Isothermal, but varies from 550® on the flanks of the rift to 300® 
within the rift, with abrupt 50 km transitions between. The transitions 
are centered on about 400 and 600 km along the section line. The model 
gives lowered heat flow within the rift, as anticipatad, which Is clearly 
In disagreement with measurements.- This Is only an Interesting exercise, 
however, because at present very little Information Is available on the 
extent of the low Curie points (laterally and vertically) and the location 
and widths of the transition zone to constrain such a model. 

Also shown In Figure 6 are measured heat flow values within 100 km of 
the section line and the heat flow profile giw^n by Figure 2, which Is 
based on averages. The measurements scatter so badly that it is difficult 
to judge the models, although predicted heat flow of model A Is probably 
too low, both for the rift and the Colorado Plateau. 

Figure 7 Is from the southern Rio Grande Rift (Section 4), where 
measured heat flow tends to be higher on the average than further north. 
This is reflected in the Am distribution and in the resulting thermal 
models. 


S'. 


I 
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3.0 PART 2; INVERSION OF CONTINUED AEROMAGNETIC DATA 

This part of the project was an investigation of whether aeromagnetic 
data could be used to Infer a Curie surface by a method analogous to that 
used with satellite data on a larger scale. The data set for this study 
was from the central Oregon Cascades flown by Oregon State University. It 
was chosen because 1) the data Is readily available on tape from the NOAA 
data center In a form suitable for easy grid averaging, 2) It Is In a geo- 
thermal 1y active region, 3) It Is of sufficient extent for valid analysis, 
and 4) a Curie depth stuoy has already been carried out on this data set by 
Connard (1979). Connard carried out depth-to-magnetic-bottom computations 
by operations in the two-dimensional frequency domain, using methodology 
made popular by a number of previous Investigators, notably Spector and 
Grant (1970). Certain authors have applied the frequency-domain approach 
to finding the configuration of the bottom of a magnetized lay<““, for 
example Gerard and Oebeglla (1975) and Hahn et a1 (1976). It seemed that 
Connard's results would provide an interesting comparison with those of the 
present study. 

3.1 Approach 

The oata was gridded, with grid spacing equal to ti. . flight line spac- 
ing (1 mile ® 1.61 km). A 76(x) by 80(y) matrix resulted. It was found 
that not all the data was present on the tape, and a block of values had to 
be filled in by hand, using Connard's original map (Figure 8). With this 
labor and elimination of some bad values a nice contour map was produced 
(Figure 9) which compares favorably with Connard's, This data was to be 
upward continued, and to reduce edge effects the grid area was extended by 
si:r.?ly folding the data over In a 25 km strip around the border (Figure 
10 ). 


Our hypothesis is that at altitudes which are "large" with respect to 
the thickness of the magnetic crust, discrete "local" sources will be mini- 
mal, that the data will predominantly reflect variations 1n magnetic crust- 
al thickness, and that to a good approximation the magnetization of the 
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layer will appear to be uniform. We imagine a Curie surfac<» level undulat- 
ing about some Intermediate level which we take as the bottom of an equiva- 
lent layer. We determine a magnetization variation In the equivalent layer 
by Inversion, then transform this to thickness variations In a layer of 
constant magnetization under the assumption that the fields due to two such 
models are virtually indlstlnqulshable In the data. 

Rather than approximating the source model with dipoles, the equiva- 
lent layer was subdivided into several arrays of vertical rectangular 
prisms, with block widths determined by the number specified In each direc- 
tion. Several Inversions were then made for the several arrays to deter- 
mine the resolution achievable, by plotting the standard deviation of the 
magnetization values In the solutions against the block spacing (or number 
of blocks in each dimension). 

Connard's results suggest shallovs Curie depths of around 6-14 km. We 
therefore selected an equivalent layer thickness of 10 km. The data was 
upward continued to 25 km (Figure 11) which experience suggests is a more 
than sufficiently "large" altitude. The continued aata was then inverted 
to a series of solutions. 

3.2 R esults 

The results of the inversions are given in Table 2. A plot of param- 
eter standard deviation against number of blocks in each direction is shown 
In Figure 12. It suggests that the unstable region is entered at block 
spacing of about 25 km (5 x 5 array of blocks). 

To convert block magnetization values in the equivalent layer to 
thickness (Curie depths), we reason as follows. The vertical integral of 
magnetization in a block (i) of variable thickness h^ and constant mag- 
netization V is uh,. The vertical integral of magnetization in the 
equivalent layer of thickness H we take to be H (u+^i), where the 
are the magnetization values solved for. We equate the two, and if H 
and u can be estimated, the can be found. The results obviously 
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depend on some Independent estimate of overall Curie depth level, either 
from spect'*al estimates or from some thermal "•^del. For u we can go two 
ways. We can assume a value and see If reasonable Curie depths result. We 
took a second approach. Connard's results Indicate a range of depths of 6 
to 14 km. We assumed that the minimum /vn In the solution was associated 
with a 6 km Curie depth, the maximum with the 14 km level, and th.<s were 
able to solve for u ; then we could ask whether It was a reasonable 
valve. For the 5x5 array u Is about 1 A/m, a reasonable value. For, 
say, the 8x8 array w ■ 6 A/m, a value commonly associated with mafic 
plutons or flows, but perhaps high for crustal rocks genera'. ly. 

The depths computed for a 5 x 5 array are shown In Figure 13. Clear- 
ly, they are In agreement with the continued field variations Figure 14 
shows the depths for the 3x3 solution for comparison with Connard's 
results, since he subdivided the area into 3x3 blocks and carried out 
spectral analysis on each separately, although only four of the nine spec- 
tra showed the spectral peak required for a depth estimate. His results 
are reproduced in Figure 15. The comparison ^s only fair. In particular, 
the continuation of the central zone of shallow depths to the north sug- 
gested by Figure 15 is not supported by Figure 14. On the other hand, both 
maps predict relatively deep magnetic bottom on the left (west) side. Both 
methodologies are baseo on simplifying assumptioiis and considerable treat- 
ment of the data; further comparisons of them will be possible when heat 
flow data from the survey area is available. 




BISINESS 4/VJ TECHNOLOCIC4L SYSTEMS, IfK, 


REFERENCES 


Ander, M. E., Geophysical study of the crust and upper mantle beneath the 
central Rio Grande Rift and adjacent Great Plains and Colorado 
Plateau. Los Alamos Scientific Laboratory Report LA-8676-T pn., 
1981. 

Bhattacharyya, 6. K., and L. K. Leu, Analysis of magnetic anomalies over 
Yellowstone National Park: Mapping of Curie point Isothermal surface 

for geothermal reconnaissance, J. Geophys. Res., 80, 4461-4465, 1975. 

Blackwell, 0. D., The thermal structure of the continental crust. In The 
Structure and Physclal Properties of the Earth's Crust, Geophys. 
Monogr. Ser., vol. 14, edited by J. 0. Heacock, pp. 169-184, AGU, 
Washington, DC, 1971. 

Blackwell, D. 0., Heat flow and energy loss in the western United States, 
Mem, Geol. Soc. Am. 152, 175-208, 19/9. 

Connard, G, G., Analysis of aeromagnetic measurements from the central 
Oregon Cascades. M. $. Thes'‘», Oregon State University, 1979. 

Cordell, L., Regional gephysical setting of the Rio Grande Rift, Geol. 
Soc. Am. Bull., 89, 1073-1090, 1978. 

Diment, W. H., T. C, Urban, J, H. Sass, B. V. Marshall, R. J. Munroe, and 
Av H. Lachenbruch, Temperatures and heat contents based on conductive 
t^'ansport cf heat, U. S. Geological Survey Circ. 726, 84-103, 1975. 

Eaton, G. P., R. R. Wahl, H. J. Prostka, 0. R. Mabey, and M. D. Kleinkopf, 
Regional gravity and tectonic pattern': Their relation Lo late 

Cenozoic epeirogeny and lateral soreadlng in the western cordillera, 
Mem. Geol. Soc. Am. 152, 51-91, 1979. 

Gallllier, S. C. and M. A. Mayhew, On the Possibility of detecting large- 
scale crustal remanent magnetization with Magsat vector magnetic 
anomaly data. Geophysical Research Letters 9, 325-328, 1982. 

Gish, 0. M., G. R. Keller, and M. L. Spar, A refraction study of deep 
crustal structure in the Basin and Range - Colorado Plateau of eastern 
Arizona, J. Geophys. Res. 86, 6029-6038, 1981. 

Hahn, A., E. G. Kind and D. C. Mishra, Depth estimation of magnetic sources 
by means of Fourier amplitude spectra. Geophysical Prospecting 24, 
237-308, 1975. 

Keller, G. R., L. W. Braile, and P. Morgan, Crustal structure, geophys- 
ical models and contemporary tectonism of the Colorado Plateau, 
Tectonophysics 16, 131-147, 1979. 



Bisi\Ess 4/<iD Technological Systems; Inc 


Lachenbruch, A, H,, and J. H. Sass, Heat flow in the United States and the 
thermal regime of the crust, in The Nature and Physial Properties of 
the Earth's Crust, Geophys. Monogr. Ser., Vol. 20, edited by J. G. 
Heacock, 625-675, AGU, Washington, O.C., 1977. 

Mayhew, M. A., Initial equivalent layer magnetization model for the United 
States based on Magsat data. EOS 62, 270, 1981. 

Mayhew, M. A., An equivalent layer magnetization model for the United 
States derived from satellite-altitude magnetic anomalies. Journal of 
Geophysical Research 87, 4837-4845, 1982a. 

Mayhew, M. A., Application of satellite magnetic anomaly data to Curie 
isotherm mapping. Journal of Geophysical Research 87, 4846-4854, 
1982b. 

Mayhew, M. A., Magsat anomaly field Inversion and interpretation for the 
U.S. Final Report, NASA Contract NAS5-26328, 1982c. 

Mayhew, M. A., B. 0. Johnson and R. A. Langel , An equivalent source model 
of the satellite-altitude magnetic anomaly field over Australia. 
Earth and Planetary Science Letters 51, 189-198, 1980. 

Mayhew, M. A. and R. H. Estes, Simulations and software development related 
to Magsat data. Final Report, National Aeronautics and Space 
Administration/Goddard Space Flight Center Contract NAS 5-25720, 
Volume 1, 87 pp., 1980. 

Mayhew, M. A. and R. H. Estes, Constrained inversion of satellite magnetic 
anomaly data. Final Report, National Aeronautics and Space 

Administration/Goddard Spice Flight Center Contract NAS 5-25698, 
Volume 3, 68 pp., 1981. 

McGetchin, R. R., and L. T. Silver, A crust-upper mantle model for the 
Colorado Plateau based on observations of crystalline rock fragments 
in the Moses rock dike, J. Geophys. Res., 77, 7022-7037, 1972. 

Pakiser, L. C. and J. N. Brune, Seismic models of the root of the Sierra 
Nevada, Science 210, 1088-1094. 

Roy, R. F., D. 0. Blackwell, and E. R. Decker, Continental heat flow, in 
The Nature of the Solid Earth, edited by E. C. Robertson, pp. 506-543, 
McGraw-Hill, New York, 1972. 

Sinno, Y. A., G. R. Keller and M. L. Spar, A cristal seismic refraction 
study in west-central Arizona. J. Geoph. Res. 86, 5023-5038, 1981. 

Shuey, R. T., D. K. Schel linger, E. H. Johnson and L. B. Alley, 
Aeromagnetics and the transition between the Colorado Plateau and the 
Basin F - ,je province. Geology, 1, 107-110, 1973. 



Business 4/<iD Technological Systems, Inc 


Shuey, R. T., 0. K. Schellinger, A. C. Tripp, and L. G, Alley, Curie depth 
determination from aeromagnetic spectra, Geophys. J. R. Astron. Soc. 
50, 75-101, 1977. 

Spector, A. and F. S. Grant, Statistical models for Interpreting aeromag- 
netic data. Geophysics 35, 293-302, 1970. 

Swanberg, C. A., and P. Morgan, Heat flow map the United States based on 
silica geothermometry. In Physical Properties of Rocks, vol. II-2 
edited by Toulouklan, H. S., W. R. Judd, and R. F. Roy, In press, 
McGraw-Hill, 1982. 

Thompson, G. A. and M. L. Zoback, Regional geophysics of the Colorado 
Plateau. Tectonophysics 35, 293-302, 1970. 

Wasllewski, P. J. and E. Padovani, Magnetization in the lower crust 
(abstract), EOS Trans. AGU 61, 222, 1980. 

Wasllewski, P. J., H. H. Thomas, and M, A. Mayhew, The Moho as a magnetic 
boundary, Geophys. Res. Lett., 6, 541-544, 1979. 



Parameters used in Curie depth/thermal 
models of Figures 4-7 
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MIN = minimum value 

MU = computed magnetization 
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FIGURE CAPTIONS 


Figure 1. Apparent magnetization contrast in 40 km thick layer derived by 
inversion of Magsat anomaly data. Contour interval 0,2 A/m, From Mayhew 
(1982c). Four section lines of Figures 4-7 shown. BR is Basin and Range, 
SR is Snake River Plain, CP is Colorado Plateau, R is Rio Grande Rift, S is 
Sierra Nevada, Wb and Ob are Williston and Denver Basins, P is western 
boundary of Great Plains. 

Figure 2. Contours of heat flow averaged over two staggered 200 km grids 
and guided by Am distribution of Figure 1. Contour Interval 0.5 hfu. 
From Mayhew (1982c). 

Figure 3. Estimates of magnetic crustal thickness based on aeromagnetic 
spectra, depth to 550®C (assumed Curie temperature) from thermal models, 
or, for Sierran region, seismic crustal thickness. 1 is Sierras along 
38®N, 2 is Sierras as a whole, 3 is Uinta Basin, 4 is Utah High Plateaus, 5 
is Yellowstone, 6 is Colorado Plateau, 7 is southwestern Great Plains, 8 is 
Basin and Range as a whole, 9 is Rio Grande Rift. Data from Roy et al 
(1972), Pakiser and Brune (1980), Eaton (1963), Lachenbruch and Sass 
(1977), Shuey et al (1977), Cordell (1978), Bhattacharyya and Leu (1975), 
McGetchin and Silver (1972), Keller et al (1979), Stewart and Pakiser 
(1962), Blackwell (1971), Cook et al (1978,79), Decker and Smithson (1975), 
and Olsen et al (1979). 

Figure 4. Section 1. Crustal seismic results from Eaton (1963), Cordell 
(1978), Jackson and Pakiser (1965), and Jackson et al (1963). Numbers are 
uppermost mantle velocities. M is Moho. Circles are heat flow values 
within about 100 km of the section (Grim et al, 1977). Heavy lines within 
crustal section labeled A,B are 550®C isotherms inferred from Am curve as 
described In text. Light lines above labeled A, 8 are associated surface 
heat flow profiles computed from thermal model. 

Figure 5. Section 2. Crustal seismic results from Warren (1969) and 

Roller (1965). Numbers and symbols as in Figure 4. 

Figure 6. Section 3. Crustal seismic results from Gish et al (1981), 

Sinno et al (1981), Olsen et al (1979), and Stewart and Pakiser (1962). 
Other numbers and symbols as in Figure 4. 

Figur. 7. Section 4. Numbers and symbols as in Figure 4. 

Figure 8. Cascades aeromagnetic survey (from Connard, 1979). 

Figure 9. Contour map of digitized Cascades data set used in present 

study. Units: Hundreds of nT, contour interval 200nT. 

Figure 10. As Figure 9, with border of extended data. 

Figure 11. Data of Figure 10 upward continued to 25 km. Units: nT, 

contour interval 5nT. 
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Figure 12, Magnetization parameter standard deviation plotted against 
block width. Numbers beside dots indicate number of blocks in each 
direction within survey area. 

Figure 13. Computed depths to magnetic bottom superimposed on Figure 11 
data for original map area. 5x5 solution. 

Figure 14. As Figure 13 for 3x3 solution. 

Figure 15. Magnetic depth results and interpretation of Connard (1979). 
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Map of the Central Oregon Cascades area showing 
areas of shallow Curie point isotherm and 
regional structural interpretation of aero- 
magnetic data. 


Flsure 15 (Connard, 1979) 





